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Abstract. In a recent work Chiu proposed to modify domain wall fermions 
that allow in an optimal way fewer number of flavors than in the standard 
case. This is done using a variant of doamin wall fermions, the so-called 
truncated overlap fermions. In this note I discuss the possibility to implement 
his proposal for the original variant of domain wall fermions. I make also 
some remarks on dynamical simulations with ultraviolet suppressed domain 
wall fermions. 
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Domain wall and overlap fermions. I will briefly reveiw t he definitions of do¬ 
main wall fermions | Kaplan, 1992 , Furman and Shamir, 1995 1. Let N be the size 
of the extra dimension, D\y the Wilson-Dirac operator, and m the bare fermion 
mass. Then, the domain wall fermion action is given by: 


( 0 . 1 ) 


N 

Sdwf ■= - 1 )^ + P+^i+i + 

2=1 


together with boundary conditions: 

P+(i>N +1 


( 0 . 2 ) 


mip i) = 0 
P-(ipo + mip N ) = 0 


Here M is the five-dimensional fermion matrix of the theory and D'' = M — Dw 
with M £ (0, 2) being a mass parameter. The alternative theo ry of chiral fermions 
on the lattice is formulted by [ Narayanan and Neuberger, 1993 j. In the domain wall 
setting this theory can by eq iuvalently fomulated using truncated overlap fermions 
|Borici, 2000a . Boriqi, 20001o| : 

(0.3) 

N 

Stov ■= ^T>tov x P = '^2'4>i[a 5 (D^~t)ilJi + ( a 5D^ + l)P + 'il)i +1 + (a 5 D^ + t)P-il; i - 1 } 

i -1 

with the same boundary conditions as above. In both cases the four dimensional 
lattice spacing is set to one. 

It is well-known that domain wall fermions in both variants can be dimen- 
sionally reduced t o give an effective theory of f our dimension al chiral fermions 
|Neuberger, 1998b , Kikukawa and Noguchi, 199S| , Boriqi, 1999a I . The Dirac oper - 


ator of the theory is given by the Neuberger overlap operator flNeuberger, 1998a| : 


(0.4) D = fl - sgn (H) 

where H = 75in case of truncated overlap fermions with transfer matrix 
|Boriqi, 2000a : 

1 + a§H 


(0.5) 


T = 


1 — (15# 


For domain wall fermions it was shown that [ Boriqi, 1999a : 

Dll „ 1 


( 0 . 6 ) 


ft = 75; 


= H- 


2-a 5 Dll 2-a 5 Dll 
Recently, Chiu proposed to use a transfer matrix which is depends on the 5D 


time slice [Chiu, 20021: 
(0.7) 


T = 

( 7, 


1 + LOidsH 


1 — 

where 07,* = 1,... ,N are positive coefficients and 05 = 1. These coefficients are 
chosen such as to correspond to the optimal Chebyshev rational approximation 
of the sign function, the Zolotarev approximation polotarev, 1877 1. In his paper, 
Chiu showed how to implement this idea in case of truncated overlap fermions. 
Below I discuss Chiu’s idea in the case of domain wall fermions. 
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Quasi-optimal domain wall fermions. I postulate the optimal domain wall 
fermion action given by: 

N 

( 0 . 8 ) - 1 )^ + P+Tpi +1 + P-ipi-i] 

i =1 

To show that such domain wall fermions are optimal it is sufficient to calculate 
tranfer matrices = 1 from M and find that they are in the optimal 

form S The fermion matrix can be written as an N x N block partition along 
the fifth dimension: 

— 1 P+ —mP- ^ 

M= P- uj 2 a 5 D^-t 

P+ 

\ — mP + P_ — 1/ 


Multiplying from the right by the permutation matrix: 

( P + P~ \ 

P+ ■■■ 

P_ 

\P- P+J 

I obtain: 

( {uJia^HPj r — 11)(P+ — mP-) uJia^HP- + 11 ^ 

U>2Cl5HP + — t 

UJN-ICL 5 HP- + 1 

+ 1)(P- — rnP+) u>Na^HP + — 11 J 

Further, multiplying this result from the left by the inverse of the diagonal matrix: 

(u]ia§HP + — 11 \ 


75 


U>2CL5HP + — 11 


I get: 

(0.9) 


T(to) := 


/ P+ — mP- —T\ 

l 


\-T n {P_ - mP+) 

with transfer matrices being defined by: 

Ti = 1- rrp (l + c+asiPP-), 

11 — Wia 5 .fi P+ 

By requiring transfer matrices being in the form: 

11 + 


u>Na,$HP + — 1 ) 

\ 


—Tn-i 

1 ) 


i = 1,..., N 


Ti = 


1 — Uid^H 
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it is easy to see that [ Borigi, 1999a ]: 

(0 - 10) H = H 2^km 

Note that TL depends in the time slice i. By ignoring the denominator the result¬ 
ing transfer matrix will be of the optimal form ( 0 )- This may lead to suboptimal 
approximation of the sign function and must be studied numerically. 

Finally to derive the four dimensional Dirac operator one has to compute: 

detDl m = ^IM 

detT(l) 

It is easy to show that the determinant of the TV x N block matrix T(to) is given 
by: 

detT( to) = det[(P+ — toP_) — T) • • • Tjy(P_ — mP + )\ 
or 

"1 777 "I — 777 

det T(m) = det[—-— 75(1 + T x • • • T N ) H--— (1 - Ti • • • T N )] 

which gives: 


D (n) = Ltm 1+ !- m 


1 + T 1 ---T N 


In the large N limit one gets the Neuberger operator (3.4). 


Ultraviolet suppressed domain wall fermions. One difficulty of dynamical 
simulations with domain wall fermions is the presence of many heavy fermion modes 
which are simulated together with the interesting light fermion. The problem, 
generally knwon as the ultraviolet slowing down becomes more severe in this case. 
A possible solution is to subtract the contribution of heavy modes in the first place 
|lzubuchi, RBC Collaboration, 2002 |. 

Here I discuss another solution which is based on the recent proposal of ultravio¬ 
let suppressed fermions [ Borigi, 2002a j. This proposal simply states that the lattice 
Dirac operator can be modified to give: 

(0.11) P = ^T 5 tanh a5r5DDWF 

a.5 /r 

where Ddwf is the input domain wall operator, a the lattice spacing and /i > 0 is a 
dimensionless parameter and ^ = 75 P 5 . P 5 is the reflection operator with respect 


to the 5-dimensional midpoint [Blum et al, 20001. The Dirac operator converges to 
the input operator in t he contimuum limit and is local and unitary as shown in 
detail in | Borigi, 2002a |. The input theory is also recovered in the limit gt —+ 00 . 
For n —* 0 one has D —> fi, i.e. a quenched theory. 

Note that the parameter p has a special role in this theory. Fixing /i at small 
values, typically fj, ~ 1 one can achieve the desired suppression of heavy fermion 
modes of the input domain wall operator. This way, one localizes the interesting 
light mode. 

At the first sight, it looks that the introduction of a matrix function can lead 
to potentially difficult computationas. In fact, it can be shown that using Lanc- 
zos based methods and multi-shift conjugate gradient iterations the computational 
overhead is small (lj. 

The benefit using the ultraviolet suppressed operator is hence twofold: one stays 
within the conventional domain wall formalism of chiral fermions and accelerates 
dynamical simulations by suppressing both heavy and ultraviolet modes. 
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End note. This work has been prepared during October 2002 when additions in 
the paper of Chiu2002 had not already appeared. It is nice to observe that his 
numerical computations support the proposal of using quasi-optimal domain wall 
fermions. 
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